Influence of uniaxial tensile stress on the mechanical and piezoelectric properties of 

short-period ferroelectric superlattice 
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Tetragonal ferroelectric/ferroelectric BaTiOs/PbTiOa superlattice under uniaxial tensile stress 
along the c axis is investigated from first principles. We show that the calculated ideal tensile 
strength is 6.85 GPa and that the superlattice under the loading of uniaxial tensile stress becomes 
soft along the nonpolar axes. We also find that the appropriately applied uniaxial tensile stress 
can significantly enhance the piezoelectricity for the superlattice, with piezoelectric coefficient 0I33 
increasing from the ground state value by a factor of about 8, reaching 678.42 pC/N. The underlying 
mechanism for the enhancement of piezoelectricity is discussed. 
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I. INTRODUCTION 

Ferroelectrics which can convert mechanical to elec- 
trical energy (and vice versa) have wide applications in 
medical imaging, telecommunication and ultrasonic de- 
vices, the physical properties of which arc sensitive to 
external conditions, such as strain, him thickness, tem- 
perature, electric and magnetic helds [H, 0, [1[ . BaTi03 
(BTO) and PbTi0 3 (PTO), as prototype ferroelectric 
materials and simple systems, have been intensively stud- 
ied 0, [j| . It is known that the ferroelectricity arises from 
the competition of short-range repulsions which favor the 
paraelectric cubic phase and Coulomb forces which fa- 
vor the ferroelectric phase @, 0|- As the pressure in- 
creases, the short-range repulsions increase faster than 
the Coulomb forces, leading to the reduced ferroelectric- 
ity. Accompanied with the suppression of ferroelectric- 
ity, the piezoelectricity decreases and even disappears. 
However, recent studies have shown that the noncollinear 
polarization rotation, occurring at phase transition pres- 
sure, can result in the giant piezoelectric response [HQ. 
In contrast to previous theoretical studies of the effects 
of epitaxial strain on the spontaneous polarization of fer- 
roelectric thin films, we have systematically studied the 
influence of uniaxial and in-plane epitaxial strains on the 
mechanical and piezoelectric prop erties of perovskite fer- 
roelectrics [13, HE [H, 0, [H, [Hf ■ So far, there has been 
no previous work on the effect of uniaxial tensile strains 
on the mechanical and piezoelectric properties of short- 
period BTO/PTO superlattices. 

Ferroelectric superlattices composed of alternating epi- 
taxial oxides ultrathin layers are currently under in- 
tensive study due to their excellent ferroelectric and 
piezoelectric properties jl6| |. Ferroelectricity can be in- 
duced in AB1O3/AB2O3 superlattice in spite of the 
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paraelectric nature of AB1O3 and AB2O3. This is 
because the coincidence of the positive and negative 
charge centers is destroyed in the superlattice and elec- 
tric dipoles are induced. Moreover, ferroelectricity can be 
enhanced in ferroelectric superlattices in certain stack- 
ing sequences fvi\ . The overall polarization of three- 
component SrTi0 3 (STO)/BTO/PTO ferroelectric su- 
perlattices can also be improved by increasing the num- 
ber of BTO and PTO layers d8|. Thanks to the periodic 
nature, it is possible to study the effect of uniaxial or 
biaxial strains on the properties of ferroelectric superlat- 
tices from first principles. 

In this work, we perform total energy as well as linear 
response calculations to study the effect of uniaxial ten- 
sile stress along the c axis on the mechanical and piezo- 
electric properties of short-period BTO/PTO superlat- 
tice. We show the mechanical properties by calculating 
the ideal tensile strength, elastic constants and valence 
charge density at different strains. We also show the in- 
fluence of uniaxial stress on the piezoelectricity. To reveal 
the underlying mechanisms, we study the effects of uni- 
axial tensile stress on the atomic displacements and Born 
effective charges, respectively. 



II. COMPUTATIONAL METHODS 

Our calculations are performed within the local den- 
sity approximation (LDA) to the density functional the- 
ory (DFT) as implemented in the plane-wave pseudopo- 
tential ABINIT package [19]. To ensure good numerical 
convergence, the plane-wave energy cutoff is set to be 
80 Ry and the Brillouin zone integration is performed 
with 6x6x6 k-meshpoints. The norm-conserving 
pseudopotentials generated by the OPIUM program are 
tested against the all-electron full-potential linearized 
augmented plane- wave method 0, [5lJ • The orbitals of 
Ba 5s 2 5p 6 6s 2 , Pb 5d 10 6s 2 6p 2 , Ti 3s 2 3p 6 3d 2 4s 2 , and O 
2s 2 2p 4 are explicitly included as valence electrons. The 
dynamical matrices and Born effective charges are com- 
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puted by the linear response theory of strain type per- 
turbations, which has been proved to be highly reliable 
for ground state properties [H, [H, HH • The polariza- 
tion is calculated by the Berry-phase approach [25[ . The 
LDA is used instead of the generalized gradient approx- 
imation (GGA) because the GGA is found to overes- 
timate both the equilibrium volume and strain for the 
perovskite structures [26| . The piezoelectric strain coeffi- 



cients cb, 



L^S/ii/, where e is the piezoelectric stress 



tensor and the elastic compliance tensor s is the recipro- 
cal of the elastic stiffness tensor c (Roman indexes from 
1 to 3, and Greek ones from 1 to 6). 

In the calculations, a double-perovskite ten-atom su- 
percell along the c axis is used for the tetragonal short- 
period BTO/PTO superlattice. The primitive periodic- 
ity of tetragonal structure with the space group P4mm 
is retained, which is more stable in energy than the 
rhombohcdral structure. For the tetragonal perovskite 
structure compounds BTO and PTO, the equilibrium 
lattice parameters are a(BTO)=3.915 A, c(BTO)=3.995 
A, a(PTO)=3.843 A, and c(PTO)=4.053 A, which are 
slightly less than the experimental values of 3.994, 4.034, 
3.904, and 4.135 A, respectively [H, G3. A sketch of 
ground-state short-period BTO/PTO superlattice with 
its atomic positions is shown in Fig.l. 




Ba {0. 
Pb (0. 
Ti1 (0. 
Ti2 (0. 

01 (0. 

02 (0. 

03 (0 

04 (0. 

05 (0 

06 (0. 

[001] 



0.0 0.000) 
0.0 0.497) 
5 0.5 0.249) 
5 0.5 0.741) 
5 0.0 0.262) 
0.5 0.262) 
5 0.5 0.012) 
5 0.0 0.755) 
0.5 0.755) 
5 0.5 0.508) 



.[010] 



.[100] 



FIG. 1: The sketch of short-period ferroelectric superlattice 
with its atomic positions. 

To calculate the uniaxial tensile stress 033, we apply a 
small strain increment 773 along the c axis and then con- 
duct structural optimization for the lattice vectors per- 
pendicular to the c axis and all the internal atomic posi- 
tions until the two components of stress tensor (i.e., an 
and 022) are smaller than 0.05 GPa. The strain is then 
increased step by step. Since <7ii=7?i(cii + C12) + 773C13, 
the elastic constants satisfy 773/771 ~-(cn +c 12 )/c 13 under 
the loading of uniaxial tensile strain applied along the c 
axis, where the strains r\i are calculated by 771 = 772 = 



(a — ao)/ao and 773 = (c — cq)/cq, with ao=3.897 A and 
co=7.859 A being the lattice constants of the unstrained 
superlattice structure. We have examined the accuracy 
of our calculations by studying the influence of different 
strains on the properties of BTO and PTO, respectively 

mm qui. 



III. RESULTS AND DISCUSSION 
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FIG. 2: (a) Uniaxial tensile stress as a function of tensile 
strain 773, and the inset reflects the relation between strains 
773 and 771. (b) Elastic constants as a function of stress CT33. 



Figure 2(a) shows the uniaxial tensile stress 033 as a 
function of strain 773 . The relation between strains 771 and 
773 is shown in the inset, which satisfy 773 > —2771. The 
stress 1T33 increases until reaching its maximum value of 
6.85 GPa with increasing strain, indicating that the cal- 
culated ideal tensile strength is 6.85 GPa for the superlat- 
tice, which is the maximum stress required to break the 
superlattice. Figure 2(b) shows the elastic constants as a 
function of stress 033, which reflect the relation between 
stress and strain. What is the most unexpected is that 
the constant C33 first decreases until reaching its mini- 
mum value at <r c =3.26 GPa and then gradually increases, 
promising a large electromechanical response at a c [27| . 
The minimum C33 corresponds to the minimum slope of 
the curve of Fig. 2(a) at a c . Other elastic constants, espe- 
cially en, always decrease with increasing 033, indicating 
that the superlattice under the loading of uniaxial tensile 
stress along the c axis becomes soft along the nonpolar 




FIG. 3: Calculated valence charge density along the c axis 
in the (100) (a) and (200) (b) planes of superlattice at equi- 
librium (1), maximum piezoelectric coefficient (2) and ideal 
tensile strength (3). 



To illustrate the change of chemical bonds with uni- 
axial tensile stress, figures 3(a) and 3(b) are plotted to 
show the valence charge density along the c axis in the 
(100) and (200) planes of the superlattice at equilib- 
rium, maximum piezoelectric coefficient and ideal ten- 
sile strength, respectively. The Pb-C>4, Ba-Oi, Tii-Oi 
and Ti2-04 bond lengths are not sensitive to the uniaxial 
tensile stress along the c axis, suggesting that the orbital 
hybridizations between these atoms are not sensitive to 
the uniaxial strain, whereas the Tii-03 and T12-O6 bonds 
elongate remarkably with increasing stress. Following the 
evolution of the charge density, we find that the weak 
Tii-03 bond starts to break first, followed by the Ti 2 -06 
bond. After the bond breaks, the system converts into 
a planar structure with alternating layers. On the other 
hand, figures 3(a) and 3(b) show that the valence charge 
density becomes more and more unsymmetrical with the 
uniaxial tensile stress increasing, indicating the increase 
in polarization. To confirm this, we have directly cal- 
culated the relations between the polarization and the 
uniaxial tensile stress with the Berry-phase approach. 

Figure 4(a) shows the polarization as a function of uni- 
axial tensile stress. For the ground-state superlattice, 
the calculated spontaneous polarization of 0.29 C/m 2 is 
less than the theoretical value of 0.81 C/m 2 of ground- 
state PTO, but slightly larger than the value of 0.28 
C/m 2 of te trag onal BTO (the other theoretical value is 



0.26 C/m 2 [28j), which supports the conclusion that the 




sharp interfaces suppress the polarization in short-period 
BTO/PTO superlattices [28]. As the stress 033 increases, 
the polarization dramatically increases with the maxi- 
mum slope appearing at a c , indicating that the ferroelec- 
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FIG. 4: Uniaxial tensile stress dependence of (a) polarization 
and (b) piezoelectric coefficients (i.e., dzi and 



trie phase becomes more and more stable with respect to 
the paraelectric phase. 

Figure 4(b) shows the variation of piezoelectric coef- 
ficients with stress 033, which are calculated by the lin- 
ear response theory. The piezoelectric coefficients all in- 
crease with increasing 1733 and reach their maximum val- 
ues at (7 C , indicating that the appropriately applied uni- 
axial tensile stress can enhance the piezoelectricity for the 
superlattice. The piezoelectric coefficient (I33 of ground- 
state superlattice is 86.36 pC/N, which is slightly less 
than the value of 103.18 pC/N of PTO, but much larger 
than the value of 36.43 pC/N of BTO. Under the loading 
of uniaxial tensile stress applied along the c axis, (I33 is 
increased from its ground state value by a factor of about 
8, reaching 678.42 pC/N for the superlattice. From pre- 
vious calculations [141 ] , we know that the uniaxial tensile 
stress can only enhance (I33 of PTO to the maximum 
value of 380.50 pC/N. The enhancement of piezoelec- 
tricity is supported by the conclusion of uniaxial tensile 
stress dependency of elastic constant C33 [see Fig. 2(b)]. 
Note that the polarization under uniaxial stress remains 
along the < 001 > direction and that the piezoelectric 
coefficients reflect the slope of polarization versus stress 
curves. The enhancement of piezoelectricity corresponds 
to the maximum slope of the curve of Fig. 4(a) at a c , it 
is the change of magnitude of polarization that leads to 
the enhancement of piezoelectricity. 

To reveal the underlying mechanisms for the abnor- 
mal piezoelectricity, we study the effects of uniaxial ten- 
sile stress on the Born effective charges and atomic dis- 
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FIG. 5: (a) Born effective charges Zz Z and (b) atomic displace- 
ments along the c axis (in c units) , relative to the centrosym- 
metric reference structure, as a function of uniaxial tensile 
stress. 

placements, respectively [see Figs. 5(a) and 5(b)]. Since 
the atomic displacements and polarization are all along 
the c axis, only charges Z* z contribute to the polariza- 
tion. The uniaxial tensile stress reduces the effective 
charges, which remain almost constant when 1733 > a c . 
The charges Z* z of Oi and O4 atoms are much close to 
their normal charges, so does the case of Ti atoms when 
033 > ct c , whereas Z* z of O3 and Oq atoms are anoma- 
lously large compared with their normal charges, sug- 
gesting the strong orbital hybridization between Tii (and 
Ti 2 ) 3d and 6 (and 3 ) 2p states [see Fig.3(b)]. Note 
that the Ba atom is fixed at (0, 0, 0) during the first- 



principles simulations. The displacements of O atoms, 
which are much larger than those of Pb and Ti atoms 
for a broad range of stress, are greatly enhanced as the 
stress (T33 increases, especially near a c , leading to the 
drastic increase in polarization. It is concluded that as 
the stress 033 increases, the atomic displacements are so 
greatly enhanced that the overall effect is the increase in 
polarization, even though the magnitudes of Z* z decrease 
with the stress increasing. 



IV. SUMMARY 

In summary, we have studied the influence of uniaxial 
tensile stress applied along the c axis on the mechanical 
and piezoelectric properties of short-period BTO/PTO 
superlattice using first-principles methods. We show that 
the calculated ideal tensile strength is 6.850 GPa and 
that the superlattice under the loading of uniaxial ten- 
sile stress becomes soft along the nonpolar axes. We 
also find that the appropriately applied uniaxial tensile 
stress can significantly enhance the piezoelectricity for 
the superlattice. Our calculated results reveal that it is 
the drastic increase of atomic displacements along the c 
axis that leads to the increase in polarization and that 
the enhancement of piezoelectricity is attributed to the 
change in the magnitude of polarization with the stress. 
Our work suggests a way of enhancing the piezoelectric 
properties of the superlattices, which would be helpful to 
enhance the performance of the piezoelectric devices. 
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